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ABSTRACT 

We present results on the star-formation rate (SFR) versus stellar mass (M*) relation (i.e., the 
“main sequence”) among star-forming galaxies at 1.37 <z< 2.61 using the MOSFIRE Deep Evolution 
Field (MOSDEF) survey. Based on a sample of 261 galaxies with Ha and H/3 spectroscopy, we have 
estimated robust dust-corrected instantaneous SFRs over a large range in M* 10® ® — 10®^ ® M©). 

We find a correlation between log(SFR(Ha)) and log(M,) with a slope of 0.65±0.08 (0.58±0.10) 
at 1.4 <2 <2.6 (2.1 < z < 2.6). We find that different assumptions for the dust correction, such as 
using the color-excess of the stellar continuum to correct the nebular lines, sample selection biases 
against red star-forming galaxies, and not accounting for Balmer absorption can yield steeper slopes 
of the log(SFR)-log(M*) relation. Our sample is immune from these biases as it is rest-frame optically 
selected, Ha and H/3 are corrected for Balmer absorption, and the Ha luminosity is dust-corrected us¬ 
ing the nebular color-excess computed from the Balmer decrement. The scatter of the log(SFR(Ha))- 
log(M») relation, after accounting for the measurement uncertainties, is 0.31 dex at 2.1 < z < 2.6, which 
is 0.05 dex larger than the scatter in log(SFR(UV))-log(M,). Based on comparisons to a simulated 
SFR-M, relation with some intrinsic scatter, we argue that in the absence of direct measurements 
of galaxy-to-galaxy variations in the attenuation/extinction curves and the IMF, one cannot use the 
difference in the scatter of the SFR(Ha)- and SFR(UV)-M* relations to constrain the stochasticity of 
star formation in high-redshift galaxies. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: 
star formation 


1. INTRODUCTION 

Star-forming galaxies show a relatively tight correla¬ 
tion between their star-format ion rate (SFR) and stel¬ 
lar mass. This relation, commonly known as the star¬ 
forming “main sequence” (Noeske et al. 2007a), has been 
intensively studied over the past decade (e.g., at z ~ 2, 
Reddy et al. 2006a; Daddi et al. 2007a; Pannella et al. 
2009; Rodighiero et al. 2011; Whitaker et al. 2012b; 
Wuyts et al. 2011b; Reddy et al. 2012b; Sawicki 2012; 
Schreiber et al. 2014; Whitaker et al. 2014b; Atek et al. 
2014; Rodighiero et al. 2014). The main sequence is gen¬ 
erally assumed to be a linear relation between log(SFR) 
and log(M*) with an associated scatter due to observa¬ 
tional uncertainties and intrinsic scatter. The slope of 
the relation tells us about star formation efficiency as 
a function of stellar mass (Whitaker et al. 2014b; Gen- 
zel et al. 2015). The intrinsic scatter of the log(SFR)- 
log(M*) relation is predicted to be due to variations in 
the gas accretion histories of different galaxies (Dutton 
et al. 2010a) and indicates the level of burstiness in star 
formation history. The common physical interpretation 
of the tightness of the SFR-M* relation is that the bulk 
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of the stellar mass in star-forming galaxies is built in 
a relatively steady process, as opposed to a rapid star- 
burst mode, as occurs in major mergers. The merger- 
driven starburst galaxies with high specific SFRs tend to 
lie above the main sequence relation (Rodighiero et al. 
2011; Atek et al. 2014) and the quiescent galaxies that 
have little ongoing star formation, identified by their red 
near-infrared colors (Williams et al. 2009; Nayyeri et al. 
2014), populate a region below the sequence (Noeske 
et al. 2007a; Wuyts et al. 2011b). 

It has been argued that the majority of galaxies spend 
most of their lifetime on the main sequence (Noeske et al. 
2007a). Therefore, studying the parameters of the main 
sequence and comparing them with the galaxy formation 
models can shed light on our understanding of processes 
that govern galaxy evolution, such as stellar and AGN 
feedback (Kannan et al. 2014; Torrey et al. 2014; Sparre 
et al. 2015), gas accretion rates (Dutton et al. 2010a), and 
gas inflows and outflows (Dave et al. 2011a). The SFR- 
M* relation can be used as a test for galaxy evolution 
simulations to see how well the simulations can reproduce 
the observed universe (e.g., Behroozi et al. 2013b). 

Determining the slope and intrinsic scatter of the 
log(SFR)-log(M*) relation can be hindered by sample se¬ 
lection effects, the diagnostic used to infer the SFR, and 
the method used to correct observed SFR for dust at¬ 
tenuation (for a summary see, Speagle et al. 2014). The 
UV-selected samples and the UV-inferred SFRs are bi¬ 
ased against massive and dusty galaxies where the bulk 
of star formation is obscured and the UV slope is de¬ 
coupled from extinction. Moreover, as SED- and UV- 
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based SFRs are dependent on the same stellar popula¬ 
tion models used to infer the stellar masses, the scat¬ 
ter in the log(SFR)-log(M*) relation may be underesti¬ 
mated when using SED or UV based SFR tracers (Reddy 
et al. 2012b). In order to overcome this issue, and more 
robustly measure SFRs for dustier galaxies, some stud¬ 
ies have used SFR(IR)-|-SFR(UV) as a proxy for bolo- 
metric SFR (e.g., Reddy et al. 2006a; Schreiber et al. 
2014; Whitaker et al. 2014b). However, estimating the 
SFR(IR) based solely on Spiteer/MIPS 24/rm imaging 
has its shortcomings and could artificially tighten the 
main sequence if the scatter in L(IR) at a given L(8^m) 
luminosity is not taken into account (e.g., Utomo et al. 
2014; Hayward et al. 2014). 

We improve upon previous high-redshift studies by uti¬ 
lizing a representative sample of galaxies with measure¬ 
ments of the Ha and H/3 emission lines. These lines of¬ 
fer the advantage of probing star formation on shorter 
timescales than the UV and IR, they are largely inde¬ 
pendent of the stellar population, and thus can be used 
to more accurately assess the intrinsic scatter in the 
log(SFR)-log(M,) relation. 

The star-formation activity was at its peak at z ^ 2 
(Reddy & Steidel 2009; Shapley 2011), making that red- 
shift a critical epoch to study the evolution of galaxies. 
Obtaining rest-frame optical spectra of large samples of 
galaxies at z ^ 2 has been challenging until recently due 
to the high terrestrial background in the near-IR and 
the lack of multi-object near-IR spectrographs on 10- 
meter-class telescopes. The newly-commissioned MOS- 
FIRE spectrograph (McLean et al. 2012) on the 10 m 
Keck I telescope has enabled us to build a large sample 
of galaxies at z ~ 2 with coverage in both the Ha and 
H/3 lines. In this study, we use a large sample of 261 
rest-frame optically selected and spectroscopically con¬ 
firmed galaxies at 1.4 z 2.6, which was obtained as part 
of the MOSEIRE Deep Evolution Field (MOSDEF) sur¬ 
vey (Kriek et al. 2015). To study the SFR-M* relation, 
we use SFRs based on the Ha luminosity that are ro¬ 
bustly corrected for dust attenuation using the Balmer 
decrement (L(Ha)/L(H/3)). We demonstrate how differ¬ 
ent SFR diagnostics and dust correction recipes can alter 
the slope and scatter of the log(SFR)-log(M,) relation, 
and thus potentially affect our conclusions regarding the 
way in which galaxies build up their stellar mass at high 
redshift. 

The outline of this paper is as follows. In Section 2, 
we present our sample properties and our measurements 
including SFRs and stellar masses. In Section 3, we de¬ 
scribe the methodology we used to derive the log(SFR)- 
log(M*) relation parameters. We compare the scatter 
of the log(SFR)-log(M*) relation for different SFR indi¬ 
cators in Section 4. In Section 5, we discuss the slope 
of the log(SFR)-log(M*) relation and the potential bi¬ 
ases arising from sample selection, dust correction, and 
Balmer absorption correction. We compare our findings 
to previous studies in Section 6, and finally, the results 
are summarized in Section 7. Throughout this paper, a 
Chabrier (2003) Initial Mass Function (IMF) is assumed 
and a cosmology with Hg = 70 km s“^ Mpc“^, Ha = 0.7, 
and flm = 0.3 is adopted. 

2. SAMPLE AND MEASUREMENTS 


2.1. The MOSDEF Survey 

The MOSDEF survey is a large multi-year survey with 
the MOSFIRE multi-object spectrometer on the Keck 
I telescope (McLean et al. 2012). The aim of the sur¬ 
vey is to obtain rest-frame optical spectra of ^ 1500 H- 
selected galaxies to study their stellar, gaseous, metal, 
dust, and black hole content. The MOSDEF observa¬ 
tions are conducted in five fields of the Cosmic Assem¬ 
bly Near-IR Deep Extragalactic Legacy Survey (CAN- 
DELS; Koekemoer et al. 2011; Grogin et al. 2011), con¬ 
sisting of AEGIS, COSMOS, GOODS-N, GOODS-S, and 
UDS, targeting three redshift bins: 1.37 < z < 1.70, 
2.09 < z < 2.61, and 2.95 < z < 3.80. Targets for 
spectroscopy are prioritized by their H-band magnitude 
(using the 3D-HST catalogs, Skelton et al. 2014) and the 
availability of spectroscopic, grism, and photometric red- 
shifts, down to H = 24.0, 24.5, and 25.0 magnitude, for 
each redshift bin respectively. 

Details of the survey strategy, observations, data re¬ 
duction, and characteristics of the full galaxy sample are 
described in Kriek et al. (2015). In this study, we use the 
data accumulated from the first two years of the survey. 

2.2. Ha Sample 

Out of the full MOSDEE sample we selected a sub¬ 
sample of 342 galaxies that have secure redshifts and 
coverage of both the Ha and H/3 lines. These criteria 
limited our sample to the low and middle redshift bins 
that correspond to 1.37 < z < 1.70 and 2.09 < z < 2.61, 
respectively. There were 49 AGN that were identi¬ 
fied and removed from the sample based on their IR 
and X-ray properties, as well as high [Nil]/Ha ratios 
([Nii]/Ha> 0.5, Coil et al. 2015). In order to have se¬ 
cure SFR measurements, we further limited our sample 
to objects where the Ha line was detected with signal- 
to-noise (S/N) > 3. This selection resulted in a sample 
of 264 objects. 

Among the Ha-detected galaxies, the H/3 fluxes of 19% 
(48 out of 264) were below 3cr. Out of these 48 H/3- 
undetected objects, 40 (83%) had a bright sky line close 
to the H/3 line, while for the detected H/3 sample the 
fraction was only 34%. The higher frequency of skyline 
contamination in the H/3-undetected sample contributes 
to the lower H/3 detection rate among these galaxies. 

To assess the relation between SFR and M, for star¬ 
forming galaxies, we proceeded to remove all those galax¬ 
ies that were deemed to be “quiescent”. In particular, 
we used a rest-frame U-V versus V-J color selection 
(Williams et al. 2009) to identify the quiescent galaxies. 
There was only one object in the Ha-H/3 detected sample 
and two in the H/3 undetected sample that were identified 
as quiescent and thus removed from the sample, leaving 
a final sample of 261 galaxies, with 215 detected in both 
lines. The redshift distribution is shown in Figure 1. 

2.3. Line Flux and Balmer Absorption Measurements 

For each object, an initial redshift was estimated based 
on the emission line with the highest S/N ratio. The line 
fluxes were then measured by fitting Gaussian functions 
to the line profiles. The Ha and [Nil]AA6548,6584 dou¬ 
blet were fit using three Gaussian functions. To estimate 
the uncertainty in the line fluxes, the spectra were per- 
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Figure 1. Redshift distribution of 261 star-forming galaxies in 
our sample with Ha detections and coverage of the H/3 line. There 
are 88 objects at 2 : ~ 1.5 and 173 at 2 ~ 2.3. 

turbed within the error spectra ^ and all emission lines 
were re-fit. The dispersion among perturbed spectra was 
taken as the uncertainty in the line fluxes (Kriek et al. 
2015). 

An important aspect in determining the Balmer emis¬ 
sion line fluxes is to correct for the underlying absorption 
that is produced in the atmosphere of primarily A-type 
stars. The Balmer absorption fluxes were determined 
from the best-fit SED models to the multi-wavelength 
photometry (Reddy et al. 2015), and were added to the 
emission line fluxes to get the total, absorption-corrected 
fluxes. For objects with Ha (H/3) detected at > Scr level, 
ratio of the average Ha (H/3) Balmer absorption flux to 
the average Ha (H/3) observed flux is ~ 2% (~ 15%). 
The uncertainty in the absorbed fluxes is ~ 2% (Reddy 
et al. 2015), which is negligible compared to the typical 
uncertainty of our emission line fluxes (~ 15%), and we 
decided not to include it in our flux uncertainty estima¬ 
tions. 

2.4. Mass and SFR Determination 

The stellar masses were determined through a min¬ 
imization SED fitting procedure assuming Bruzual & 
Chariot (2003) models for solar metallicity, a Chabrier 
IMF, the Calzetti et al. (2000) dust attenuation curve 
with 0.0 < E{B — V) < 0.6, rising star-format ion histo¬ 
ries® and a 50Myr lower limit on the age of the galax¬ 
ies, assumed to be the typical dynamical time scale at 
this redshift (see, Reddy et al. 2015). The SED-fitting 
was performed using photometry from the 3D-HST cat¬ 
alogs (Brammer et al. 2012; Skelton et al. 2014) that has 
been corrected for the strongest emission lines measured 
in the MOSDEF survey, including [Oii]AA3727,3730, 

^ Measuring the error spectra of the MOSDEF observations is 
described in detail in Kriek et al. (2015). In brief, we constructed 
noise spectra by three independent methods: 1. based on the sky 
and readout noise level of the detector, 2. based on the variations 
between the science frames, and 3. based on the variations between 
the reduced spectra. The measured noise spectra based on the 
three methods agree with each other very well, showing that the 
noise measurements are robust. We used the first method error 
spectra in our analysis. 

® We should note that assuming a constant star-formation his¬ 
tory for the SED modeling does not change the parameterization of 
the SFR-M* relation, as the scatter in the ratio of stellar masses in¬ 
ferred from a constant star-formation history modeling to the stel¬ 
lar masses derived assuming an exponentially rising star-formation 
history is ~ 0.06, which is smaller than the typical relative error of 
the stellar masses (~ 0.10, see Section 3.2). 


[Oiii]AA4960,5008, H^, and Ha (Reddy et al. 2015). 
Errors in the SED parameters were estimated through 
Monte Carlo simulations; each photometric flux was ran¬ 
domly perturbed assuming a Gaussian distribution cen¬ 
tered at the measured value with a standard deviation 
equal to the photometric error. This procedure was re¬ 
peated many times and each photometric realization was 
refit to find the best-fit SED model. The standard devi¬ 
ation of the simulated SED parameters was taken as the 
error in the respective parameter. 

We converted the Ha luminosity to SFR using the 
Kennicutt (1998) relation modified for a Chabrier (2003) 
IMF, and corrected for dust attenuation using the neb¬ 
ular E{B — V) derived from the Balmer decrement 
(Ha/H/3), assuming a Galactic extinction curve taken 
from Cardelli, Clayton, & Mathis (1989, CCM). To de¬ 
termine the standard deviations in the corrected SFRs, 
we perturbed the measured Ha and Bj5 fluxes according 
to their uncertainties by assuming a Gaussian distribu¬ 
tion centered at the measured value with a standard de¬ 
viation equal to the measurement error and recalculated 
the corrected SFRs based on each flux realization. We 
also accounted for spectroscopic slit losses® by adding 
an additional uncertainty of 16% (20%) to the Ha (H,5) 
measurement uncertainties in quadrature (Kriek et al. 
2015). Where H/3 is undetected, a 3cr upper limit on the 
H/3 flux was used to estimate a lower limit on the nebular 
E{B — V) and the corresponding dust-corrected SFR. 

We used the ancillary 3D-HST broad- and medium- 
band photometry at rest-wavelengths 1268 - 2580 A to 
find the UV luminosity at 1600 A and the UV slope 
by fitting a power law to the photometry (/> oc A“). 
The UV slope was used to correct the UV luminosity 
for dust attenuation, following the Meurer et al. (1999) 
relation. We did not use the SED-inferred E{B — V) 
because it is correlated with the stellar mass as both 
quantities are determined through SED fitting The 
luminosities were converted to SFRs using the Kennicutt 
(1998) relation, adjusted for a Ghabrier IMF. 

2.5. Stacking Procedure 

To further investigate the SFR-M* relation and the 
role of the H/3-undetected objects, we employed spec¬ 
tral stacks constructed in the following manner. Indi¬ 
vidual spectra were converted from flux density to lu¬ 
minosity density, normalized by the Ha luminosity, and 
interpolated to a rest-frame grid with wavelength spac¬ 
ing of 0.5 A. The composite spectrum was then created 
by taking the average of the individual spectra, weighted 
by inverse variance of the error spectra. The composite 
error spectrum was estimated as the square root of the 
variance of the weighted mean. 

® The slit loss uncertainty accounts for the loss of flux of a re¬ 
solved source outside of the slit aperture. We calculated the slit loss 
uncertainty in our flux calibration by comparing the SED-inferred 
fluxes of the continuum detected galaxies in MOSDEF with their 
spectroscopic fluxes. For a full description, see Section 3.3, Kriek 
et al. (2015). 

As it can be seen in Table 1, the SFR(SED) and the SFR(UV) 
that is dust corrected by the E{B — V) inferred from the SED mod¬ 
eling show tighter and steeper SFR-Af* relations compared to that 
of SFR(Ho) and SFR(UV) corrected by the UV slope, indicating 
a stronger correlation with M*. We should note that SFR(UV) 
corrected with the UV slope is still not completely uncorrelated 
with M* as the UV slope translates to an E{B — V). 
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Figure 2. SFR(Hq) as a function of stellar mass for star-forming 
galaxies at z = 1.37 - 2.61, known as the star-forming main se¬ 
quence. The individual points show galaxies with detected Ha and 
H/3 and the triangles are those with undetected for which 3cr 
lower limits of the dust corrected SFRs are plotted. The dashed 
line is the best linear fit to the individual detected galaxies with 
log(M*/ Mq) > 9.5. The larger circles are stacks in five mass bins, 
taking into account only the detected galaxies (the light blue cir¬ 
cles) and including the H/3 non-detections (the red circles). The 
shaded regions and numbers on the bottom represent the fraction 
of galaxies with robust redshifts and Ha detection with respect to 
the targeted objects in each mass bin. The 3cr limits on the SFR 
at z = 1.37 and 2 = 2.09 are calculated using the line sensitivities 
in H and K bands, respectively (Kriek et al. 2015), and are shown 
by dashed lines. 

The composite Ha and H/3 lines were fit with a Gaus¬ 
sian profile. The errors on the line fluxes are derived 
using the standard deviation of the distribution of line 
fluxes by perturbing the stacked spectrum according to 
the composite error spectrum and recalculating the line 
fluxes 1000 times. The composite line fluxes were cor¬ 
rected for Balmer absorption using the mean absorption, 
inversely weighted by flux error of the individual galaxies 
contributing to the bin (the individual Balmer absorp¬ 
tions were inferred from the best-fit stellar models). The 
normalized composite spectrum was used to derive the 
Ha/H/3 ratio. To assess the total flux of the lines, we 
stacked the individual spectra without normalizing them 
by the individual Ha luminosities. 

3. ANALYSIS 

The correlation between stellar mass and SFR has been 
studied intensively since Brinchmann et al. (2004). As 
a first step in studying this correlation in our Ha- and 
H/3-detected sample, we calculated the Spearmans rank 
correlation coefficient (p) between the log(SFR(Ha)) and 
log(M*). We obtained p = 0.57 at Str, indicating a sig¬ 
nificant positive correlation between the two quantities. 

The SFRs versus stellar masses of individual galaxies 
are plotted in Figure 2. The stacked SFRs in five stel¬ 
lar mass bins (the first bin is log(M*/M 0 ) < 9.5 and 
the rest are divided by 0.5 dex widths. Section 2.5) are 
also plotted on top of the individual points. In this sec¬ 
tion, we describe our low-mass completeness, the fitting 
methodology, and discuss the effect of H/3 nondetections 




Figure 3. Top: Distribution of Ft 60 IF to FS50I/P flux ratio (left) 
and the corresponding uncertainty (right) of UV-selected galax¬ 
ies in the MOSDEF parent sample with spectroscopic redshifts at 
2 = 2.09 — 2.61 (Reddy et al. 2006b) in the GOODS-N field. The 
fluxes and uncertainties are from the 3D-HST photometric catalog 
(Skelton et al. 2014). The orange histogram denotes galaxies with 
photometric redshifts that fall outside of the spectroscopic redshift 
range (bad photo-z). The blue histogram shows galaxies with pho¬ 
tometric redshifts in the same range as the spectroscopic redshift 
(good photo-z). The galaxies with bad photo-z’s have on aver¬ 
age bluer colors and larger uncertainties, indicative of a weaker 
Balmer/4000 A break compared to the good photo-z galaxies. 
Bottom: stellar mass distribution of the same samples described 
above. The orange/blue numbers indicate fraction of galaxies with 
bad/good photo-z’s in each mass bin relative to the total number of 
galaxies in the corresponding bin. The histograms in these figures 
show the fractions, not the actual numbers. Therefore, the relative 
bin heights do not correspond to the percentage numbers on the 
bottom and should not be confused with each other. 

and AGN on the SFR-M* relation. 

3.1. Mass Completeness 

As it is evident from Figure 2, the stacked SFR of the 
first mass bin is elevated compared to the best-fit line. In 
this section, we discuss the completeness of our sample 
and address the possible causes of incompleteness at low 
stellar masses. 

The redshifts of targets in the MOSDEF parent sample 
were adopted primarily from the 3D-HST grism measure¬ 
ments or other available spectroscopic campaigns (see, 
Kriek et al. 2015). When a spectroscopic redshift was 
not available, which was the case for 67% of the whole 
sample and 73% of the objects with M, < 10®’® Mq, 
redshifts were estimated using the 3D-HST photometry 
and the EAZY SED-fitting code (Brammer et al. 2008). 
Objects with secure photometric redshifts may be biased 
against those objects with weak Balmer/4000 A breaks 
(and hence younger ages and lower masses). To assess 
the magnitude of this bias, we cross correlated our sam¬ 
ple with the Reddy et al. (2006b) spectroscopic sample of 
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UV-selected galaxies in the GOODS-N field. Rest-frame 
UV-selected samples are immune to the bias mentioned 
above as their selection criteria are not dependent on the 
strength of the Balmer or 4000 A break. 

There were 81 galaxies at Zspec = 2.09 — 2.61 that were 
in common between the MOSDEF and the UV-selected 
sample. We divided the common galaxies into two sets: 
those that have photometric redshifts in the same range 
as their spectroscopic redshifts of z = 2.09 — 2.61 (good 
photo-z, 61 objects), and those with photometric red¬ 
shifts outside of this range (bad photo-z, 20 objects). 
The bad photo-z objects are not scattered symmetrically 
around Zspec = 2.09 — 2.61; their photometric redshifts 
are systematically lower than our target redshift range, 
i.e., Zphot ^ 2.09. 

The two filters that bracket the Balmer break are 
HST/WFG3 F160Wand HST/ACS F850LP, which were 
available for all 81 galaxies examined here. As expected, 
the galaxies with incorrect photometric redshifts had on 
average bluer F850LP - F160W colors, with larger un¬ 
certainties. The distributions of the flux ratios and the 
associated uncertainties are plotted in Figure 3. Were it 
not for prior spectroscopic redshifts, those galaxies with 
bad photo-z’s would have been missing from our sam¬ 
ple. The conclusion is further supported by the two-sided 
Kolmogorov-Smirnov test, which indicates a low prob¬ 
ability of P ~ 0.017 that the F850LP - F160W color 
distribution of the good and bad photo-z sets are drawn 
from the same parent distribution. 

The stellar-mass distributions of the two sets of galax¬ 
ies are shown in the bottom panel of Figure 3, with 
the fraction of galaxies in each mass bin written on the 
bottom of the plot. The distributions indicate that the 
majority of galaxies with stellar masses below 10 ®'^ Mq 
would have been missed from the MOSDFF sample if we 
had relied on their photometric redshifts. Therefore, to 
avoid any biases we limited our analysis to those with 
M* > 10® ® Mq. There are 128 and 185 galaxies above 
10 ®-® Mq at z = 2.09 — 2.61 and z = 1.37 — 2.61, respec¬ 
tively. The dashed line in Figure 2 is the linear regression 
fit to the individual galaxies with M* > 10® ® Mq. The 
low mass limit we adopted is also the limit above which 
the inferred log(SFR)-log(M*) slope (see Section 3.2) de¬ 
viates by greater than Icr from the slope over the full 
stellar mass range, including the low mass galaxies. 

The other possible bias at the low-mass end of our 
sample is a Malmquist bias, due to the fact that our re¬ 
quirement of an Ha detection may bias against objects 
with faint Ha emission. To test this possibility, we calcu¬ 
lated the fraction of galaxies that were detected in Ha to 
the number of galaxies targeted for spectroscopy, in bins 
of stellar mass. The fractions are listed on the bottom 
of Figure 2. The detection fraction of the first bin (the 
lowest mass bin) is similar to that of the higher mass 
bins, indicating that by limiting our sample to the Ha- 
detected objects we did not miss a significantly higher 
fraction of galaxies at low masses. 

The rest-frame optically selected sample of the MOS¬ 
DEF survey is less prone to missing highly obscured 
galaxies at high stellar masses, compared to the UV- 
selected samples. The effect of this potential bias of UV- 
selected samples versus our sample will be discussed in 
Section 5.2. 


3.2. Methodology 

There are various linear regression methods that can 
be used to characterize a correlation between two vari¬ 
ables. In most studies of the SFR-M* relation, the corre¬ 
lation between SFR and M* has been determined using 
the ordinary least-squares regression (OLS). In the OLS 
method, one variable (here, the SFR) is assumed as the 
dependent variable and the regression line is defined to 
be that which minimizes the sum of the squares of the 
vertical distances between the data points and the re¬ 
gression line. The OLS method should be used in science 
cases where there is an independent variable that is esti¬ 
mated confidently and will be used to predict the other 
variable. As over a galaxy’s lifetime the stellar mass is 
built relatively smoothly as opposed to the SFR, which 
can be stochastic and vary to a greater extent, we are 
more interested in determining the scatter and mean of 
SFR as a function of stellar mass. For this purpose, the 
OLS method is a valid statistical procedure to be used 
for characterizing the SFR-M* relation. 

We determined the best estimates for the OLS regres¬ 
sion slope and intercept by perturbing the data accord¬ 
ing to the stellar mass and SFR uncertainties 1000 times. 
The median of the slope and intercept distributions are 
taken as the best-fit parameters. The slope and intercept 
uncertainties were calculated based on the OLS method 
error estimates in Isobe et al. (1990). Isobe et al. (1990) 
estimated the errors of regression parameters for linear 
relations where there is an intrinsic scatter. The average 
of the uncertainties obtained for the individual realiza¬ 
tions was adopted as the uncertainty in the regression co¬ 
efficients. We found a slope of 0.65 ± 0.08 for SFR(Ha) at 
1.37 < z < 2.61, and 0.58 ± 0.10 for the 2.09 < z < 2.61 
sample only, with a lower mass limit of 10® ® Mq. The re¬ 
lation is presented in Figure 2 and Table 1. The SFR-M* 
relation was not investigated separately for the sample at 
1.37 < z < 1.70 due to the small number of galaxies at 
this redshift range. 

In the OLS method, the scatter is defined as the scat¬ 
ter in SFR at a certain stellar mass. We estimated 
the total observed scatter by calculating the 3-cr clipped 
dispersion in the distribution of vertical distances from 
the best-fit line. We propagated the mean measure¬ 
ment errors of mass and SFR^^ and subtracted them 
(in quadrature) from the total observed scatter to ob¬ 
tain the measurement-subtracted scatter. We found a 
measurement-subtracted scatter of 0.31 dex for SFR(Ha) 
at 2.1 < z < 2.6. The measurement-subtracted scatter 
still includes galaxy-to-galaxy variations in the assumed 
dust attenuation curve (Section 4). 

An alternative statistical procedure to the OLS method 
is a regression method that treats the two variables sym¬ 
metrically, such as the OLS bisector method. The OLS 
bisector is recommended to be applied in cases where an 
underlying functional relation between two independent 
variables is studied (Isobe et al. 1990). The scatter in 
this method is the scatter about the best-fit line to the 
relation. For comparison, we reported the OLS bisector 
values in Table 2. We should emphasize that as we are 
most interested in determining the scatter and mean of 
SFR at a given stellar mass, the OLS method is the one 

The uncertainty due to spectroscopic slit losses is included in 
the SFR error estimates as well (Section 2.4). 
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Table 1 

Parameters of the log(SFR) — log(M*) Linear Fit ^ 


Redshift Range 

SFR Indicator 

Slope 

Intercept Observed Scatter Measurement-Subtracted Scatter 

2 = 1.37-2.61 

Hq'^ 

0.65 

± 

0.08 

-5.40 

± 

0.86 

0.40 

0.36 

(N = 185)‘> 


0.62 

± 

0.08 

-5.03 

± 

0.80 

0.34 

0.30 


SED® 

0.80 

± 

0.05 

-6.79 

± 

0.55 

0.28 

0.25 


UVsED*^ 

0.79 

± 

0.07 

-6.58 

± 

0.67 

0.27 

0.20 

2 = 2.09 - 2.61 

Hq® 

0.58 

± 

0.10 

-4.65 

± 

1.05 

0.36 

0.31 

II 

to 

00 

O' 

UV^jd 

0.71 

± 

0.09 

-5.84 

± 

0.92 

0.31 

0.25 


SED® 

0.83 

± 

0.07 

-7.02 

± 

0.73 

0.29 

0.25 


UVsED^ 

0.75 

± 

0.09 

-6.11 

± 

0.95 

0.28 

0.18 


^ The parameters are calculated for galaxies with M* > 10®'® Mq, using an OLS linear regression method 
® Number of objects 

The Ha SFR is dust corrected by Balmer decrement, assuming the Cardelli Galactic extinction curve 
The UV SFR is dust corrected by the UV slope, assuming the Calzetti attenuation curve 
® The SFR inferred from the stellar population model, assuming the Calzetti curve for dust attenuation 
^ The UV SFR is dust corrected by the E(B — V) inferred from the SED model, assuming the Calzetti attenuation curve 


on "which we based our main conclusions. 

3.3. Nondetections 

It is important to investigate whether excluding the H/3 
nondetections from the main analysis causes any biases. 
In other words, do the H/3-undetected galaxies represent 
different physical characteristics or are they undetected 
due to observational shortcomings? We stacked the spec¬ 
tra of the H/3-undetected objects as described in Sec¬ 
tion 2.5, and derived the dust-corrected SFR(Ha) using 
the stacked spectrum. In the composite spectrum, the 
stacked H/3 line is detected at 3.9cr. The stacked SFR 
and median stellar mass of the undetected-H/3 objects 
are plotted on top of the SFR-M* relation in Figure 4 as 
the red star. The star lies below the linear fit by 2.5 ct, 
but is still within the SFR(Hq!) scatter of 0.36 dex, im¬ 
plying that these objects are not on average intrinsically 
different from the detected objects and excluding them 
should not systematically bias the SFR-M* relation. As 
mentioned in Section 2.2, this conclusion is further sup¬ 
ported by the fact that the spectra of a large fraction of 
these galaxies are contaminated by skylines. 

3.4. AGN 

As mentioned before in Section 2, AGN were identified 
through X-ray, IR, and optical line observations and re¬ 
moved from the sample. The AGN mostly reside on the 
high mass end of the main sequence and it is interesting 
to examine how they would change the slope and scat¬ 
ter of the log(SFR)-log(M*) relation. Here, we assumed 
that all of the rest-frame UV continuum emission in AGN 
host galaxies is coming from the stellar population; we 
determined the SFRs from the UV luminosity at 1600 
A, dust corrected using the UV slope, and estimated the 
masses through the SED modeling. Including the AGN 
does not change the slope significantly (0.66 ± 0.09) but 
results in a larger scatter (0.41 dex) that is mainly oc¬ 
curring at the massive end, as most of the AGN have 
larger stellar masses than the typical stellar mass of the 
star-forming objects. This is a well-understood selection 
effect for AGN that they are easier to detect in more 
massive galaxies (Aird et al. 2012). 

4. THE SCATTER OF THE LOG(SFR)-LOG(M,) RELATION 


In this section, we investigate the intrinsic scatter in 
the log(SFR)-log(M*) relation by using different SFR in¬ 
dicators and considering variations in the dust attenua¬ 
tion curve. The observed and measurement-subtracted 
scatter in SFR(Hq;) versus M* and SFR(UV) versus M*, 
calculated as described in Section 3.2, are presented in 
Table 1 and Figure 4. For a fair comparison, we limited 
our UV sample to the same sample as the Ha detected 
galaxies. 

The measurement-subtracted scatter reported here is 
affected by uncertainties such as galaxy-to-galaxy vari¬ 
ations in IMF and dust attenuation curve. In order to 
test the effect of dust attenuation variations, we simu¬ 
lated an SFR-M* relation with 1000 galaxies and an in¬ 
trinsic scatter of 0.2 dex, assuming the log(SFR)-log(M*) 
slope and intercept derived based on our MOSDEF data 
(Section 3.2). We attenuated the intrinsic SFRs to de¬ 
rive the observed (unobscured) SFRs following the re¬ 
lation, log(SFRobs)=log(SFRint) - 0AkxE{B - V). In 
this equation, the E{B — V) was calculated using a lin¬ 
ear correlation between the nebular E{B — V) and the 
intrinsic SFR with a scatter of 0.5 dex, based on the 
MOSDEF data. The reddening value, k\, was selected 
randomly according to a valid range of the most com¬ 
monly used attenuation curves in the literature; HHa 
varied between the MOSDEF (Reddy et al. 2015) and 
the Galzetti et al. (2000) curves’ reddening values, and 
«:i600 varied between the Milky Way and the SMC (Gor¬ 
don et al. 2003) reddening values. The dust-corrected 
SFRs were then reproduced by assuming a single k\ 
(the Calzetti value for the SFR(UV) and the Cardelli 
value for the SFR(Ha)). The scatter in the new dust- 
corrected log(SFR(Ha))-log(M*) relation was ~ 0.02 dex 
larger than the intrinsic simulated scatter, while the 
dust-corrected log(SFR(UV))-log(M*) scatter increased 
by ~ 0.10 dex. The larger increase in the SFR(UV) scat¬ 
ter is followed by larger variations of the attenuation 
curves at shorter UV wavelengths compared to the varia¬ 
tions of the attenuation curves at optical wavelengths^^. 
This simple simulation indicates that without knowing 

Here, we assumed that the scatter in kx does not change 
as a function of luminosity or stellar mass, but we know that the 
scatter in Kieoo larger at higher UV luminosities (e.g., Reddy 
et al. 2010, 2015). This effect would increase the recovered scatter 
in log(SFR)-log(M*) relation even to a greater degree. 
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Table 2 

Parameters of the log(SFR) — log(M*) Linear Fit, Using The OLS Bisector Method 


Redshift Range 

SFR Indicator 

Slope 

Intercept 

Observed Scatter 

Measurement-Subtracted Scatter 

2 = 1.37-2.61 

Ha= 

1.26 ± 0.07 

-11.57 ± 0.69 

0.33 

0.30 

(N = 185)® 

UV;3® 

1.14 ± 0.06 

-10.33 ± 0.62 

0.29 

0.26 


^ The parameters are calculated for galaxies with M, > 10®'® Mq 
^ Number of objects with > 10®'® Mq 

The Ho SFR is dust corrected by Balmer decrement, assuming the Cardelli Galactic extinction curve 
The UV SFR is dust corrected by the UV slope, assuming the Calzetti attenuation curve 


10^ 


10 ^ 


'a 10" 

X 

cc 

Li_ 

1/1 


= 0.36 dex 

"yr H/3 undetected, Ha detected 


10 “ 



1.37 < z < 2.61 


10 " 

Mass [Mg] 



Mass [Mg] 


Figure 4. SFR as a function of stellar mass for star-forming galaxies at z = 1.37 — 2.61. Left: SFR(Ho) - corrected for dust attenuation 
assuming the Balmer decrement and a Cardelli et al. (1989) extinction curve - versus M,. The red star shows the weighted-average stack 
of the H/3-undetected galaxies, where its mass is the median of the mass of galaxies in the stack. The stacked point is within the SFR(Hq) 
scatter of 0.36 dex, suggesting that those galaxies undetected in H/3 have on average a similar proportion of SFR given their stellar mass, 
as those galaxies with both lines detected. Right: SFR(UV) versus M,. The SFR is inferred from the UV luminosity at 1600A and is 
dust corrected by the UV slope, both determined from the multi-band photometry at 1268-2580A (Skelton et al. 2014). The green lines 
show the regression lines fitted to log(SFR) — log(M*) relation. The measurement-subtracted scatter in SFR (i.e., the scatter after the 
subtraction of SFR and mass measurement uncertainties) is provided on the upper left corners of each plot. 


the attenuation curve variations, it is not trivial to in¬ 
terpret the difference in the SFR(Hq;) and the SFR(UV) 
scatter solely as the SFR stochasticity due to different 
timescales. Variations in the assumed IMF and metallic- 
ity would also affect the SFR(Ha) and SFR(UV) scatter 
to different degrees, again making it difficult to conclude 
about the burstiness of the star formation based on the 
difference in the scatter of SFR(Ha) and SFR(UV). It is 
worth pointing out that picking a random star-formation 
history between a constant and an exponentially rising 
star-formation history affects the SED inferred stellar 
masses by < 0.04 dex, which is negligible compared to 
the scatter in the log(SFR)-log(M*) relation. 

The calculated measurement-subtracted scatter in the 
log(SFR(Ha))-log(M*) relation is 0.06 dex larger than 
that of the log(SFR(UV))-log(M*) relation. As the 
two SFR indicators trace SFR on different timescales^^, 
the difference in SFR scatter may be attributed to the 
stochasticity of SFR. It has been suggested in the galaxy 
formation simulations that for a massive galaxy (z = 0 
halo mass of 10^^ hi©) with a typical merger history, the 
dispersion in SFR averaged over timescales of ^ 10® years 
could be ^ 0.03 — 0.1 dex smaller than the scatter of SFR 

The Ho luminosity traces SFR on shorter timescales compared 
to the UV luminosity, as Ha is sensitive to the most massive (M* > 
15 Mq) and short-lived (ages < 10 Myr) stars, while UV traces stars 
with M* > 5Mq and ages < 100 Myr, (Kennicutt 1998; Kennicutt 
&; Evans 2012; Madau Dickinson 2014). 


averaged over shorter timescales of ^ 10^ years (Hopkins 
et al. 2014; Sparre et al. 2015; Dominguez et al. 2015). 
On the other hand, in our study the Ha light is cap¬ 
tured through a slit while the UV light is obtained from 
imaging data, which may introduce additional scatter to 
our measurements. Considering these uncertainties and 
poorly constrained effects of the variations in dust atten¬ 
uation curve and IMF, as discussed earlier in this section, 
we conclude that the 0.06 dex difference in scatter is not 
significant and can not be attributed solely to the SFR 
stochasticity in these galaxies. 

According to Speagle et al. (2014), the intercept of the 
log(SFR)-log(M,) relation changes as 0.12 ± 0.04 x t, 
where t is the age of the universe in Gyr. This change 
corresponds to an offset of 0.09 dex from 2; = 2.0 to 2.6. 
The offset has a negligible effect on the 0.31 dex scatter 
of log(SFR(Ha))-log(M*) relation derived here at 2.0 < 
z < 2.6, and subtracting it in quadrature will reduce the 
scatter to 0.30 dex. 

Prior to the use of multi-wavelength estimates of the 
SFR (e.g., SFR(IR)-|-SFR(UV)), it was common to as¬ 
sess the SFR-M* relation using SED-inferred SERs (e.g., 
Magdis et al. 2010; Sawicki 2012). To investigate how 
SFR(SED) would affect the relation, we also calculated 
the scatter using SED-inferred SFR and UV SFR, where 
the latter was corrected by the SED-derived E{B — V) 
rather than the UV slope. The scatter in SER is reduced 
to 0.25 and 0.20dex for SER(SED) and SER(UV), re- 



























spectively (the last two rows in Table 1). These values 
clearly underestimate the true scatter in the log(SFR)- 
log(M*) relation, because SFR(SED) and SED-inferred 
E{B — V) are highly correlated with the SED-derived 
stellar masses. 

Galaxy formation simulations suggest various drivers 
for the scatter in log(SER)-log(M») relation. Dutton 
et al. (2010a) argued that the source of the small scat¬ 
ter they derived (0.12dex) is the variance in mass accre¬ 
tion histories for halos of a given virial mass, but Dave 
et al. (2012) suggested that the more complex relation 
between the gas inflow rate and the dilution time, which 
shows how quickly galaxies can return to equilibrium af¬ 
ter being perturbed away from the main sequence also 
influences the scatter. Our results show that the larger 
scatter we found (~ 0.30 — 0.40 dex) is not only due 
to the measurement uncertainties, which indicates that 
the gas accretion and star-format ion histories are not as 
smooth as predicted by the some simulations. Also, by 
assuming an average attenuation curve and a single IMF 
and metallicity for all the galaxies contributing to the 
SFR-M* relation, the measurement-subtracted scatter is 
a lower limit on the intrinsic scatter in SFRs at a given 
stellar mass. 

5. THE SLOPE OF THE LOG(SFR)-LOG(M*) RELATION 

In this section, we focus on the slope of the log(SFR)- 
log(M*) relation for the SFR(Ha) sample. The inferred 
slope (see Section 3.2 for the methodology) is 0.58 ± 0.10 
for the z ~ 2.3 sample and when we include the z ^ 1.5 
sample the slope is 0.65 ± 0.08, at M* > 10®'® Mq (refer 
to Section 3.1 for the discussion on the low-mass limit). 
The slope of the full sample (i.e., M* > 10®'® Mq) is 
0.56 ± 0.09. 

There are many studies in the literature that param¬ 
eterize the SFR-M* relation and their reported slopes 
span in a wide range of ~ 0.3 — 1.0 ( e.g., Daddi et al. 
2007a; Dunne et al. 2009; Pannella et al. 2009; Santini 
et al. 2009; Rodighiero et al. 2011; Karim et al. 2011; 
Zahid et al. 2012; Reddy et al. 2012b; Whitaker et al. 
2012b; Sobral et al. 2014; Atek et al. 2014; Rodighiero 
et al. 2014). Various effects could lead to discrepancies 
in the slope estimates. In this section, we try to address 
three of the most significant effects by investigating them 
in our sample. 

5.1. Dust Correction 

The method used for dust correction affects the derived 
slope of the log(SFR)-log(M*) relation. In high-redshift 
studies, a commonly used dust correction method for the 
nebular lines is to multiply the E{B — V) derived for the 
stellar continuum by 2.27 factor and assume the Calzetti 
et al. (2000) attenuation curve. It has been shown in 
some previous studies that this recipe overestimates the 
corrected Ha SFRs for UV-selected galaxies (e.g., Reddy 
et al. 2010; Steidel et al. 2014; Shivaei et al. 2015). In 
the top panel of Figure 5 we compared the Ha SFRs cor¬ 
rected by E{B — V)nebuiar, Computed from the Balmer 
decrement, with those corrected by 2.27 x E{B — V)sed- 
We emphasize that this practice was done solely for 
the purpose of demonstrating how the commonly used 
dust correction recipes could potentially affect the pa¬ 
rameterization of the SFR-M* relation in the absence of 
Balmer decrement measurements. The SFRs corrected 




Figure 5. The SFR(Ho)-Mt relation, where the SFR is corrected 
for dust attenuation by different methods and attenuation curves. 
Top: The blue dots show SFR(Ha) dust corrected by the nebu¬ 
lar E{B — V) (derived from the observed Balmer decrement) and 
assuming a Calzetti curve. The red dots represent the SFR(Hq:) 
corrected by the SED inferred E{B — V) multiplied by 2.27, again 
assuming a Calzetti curve. The latter is a method that is com¬ 
monly used in high-redshift studies to correct the Ha luminosi¬ 
ties without having access to the Balmer decrement. This method 
overestimates the corrected SFRs and results in a steeper slope 
(0.90 di 0.07) compared to the nebular corrected SFRs (slope = 
0.71 di 0.09). The dashed lines show the linear fits to galaxies with 
log(M*/M 0 ) > 9.5. The yellow line is the linear fit to galaxies 
corrected by nebular E{B — V) but assuming a Cardelli Galactic 
extinction curve, which we used as a default in this study. Bottom: 
The SFR-M* relation for the same set of galaxies, corrected with 
the nebular E(B — V) assuming three different dust curves: the 
Calzetti curve, the MOSDEF curve, and the Cardelli curve. The 
best-fit lines are shown with the relevant colors. All three slopes 
are consistent with each other within their uncertainties. 

with 2.27 X E{B — V)sed assuming the Calzetti curve 
yield a main sequence slope of 0.90 ± 0.05, while the 
slope of those corrected by the E{B — V)nebuiar with 
the Calzetti curve is 0.71 ± 0.09. The steeper slope is 
expected as the color-excess {E{B — V)) increases with 
stellar mass (e.g.. Cam & Best 2010) and multiplying 
E{B — V) by a factor (here, 2.27) causes a more signif¬ 
icant difference at higher masses. Adopting a value be¬ 
tween 1 and 2.27 (Kashino et al. 2013; Wuyts et al. 2013) 
will make the slope shallower than 0.90 but still steeper 
than the 0.71 slope derived by the E[B — V)nebuiar with 
the Calzetti curve. 

We also investigated the potential effect of the assumed 
attenuation curve on the best-fit slope of the log(SFR)- 
log(M*) relation. We corrected the observed Ha SFRs 
assuming different attenuation curves: the Calzetti curve 
(Calzetti et al. 1994), the attenuation curve derived from 
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Figure 6. Top: distribution of our Ha and H/3 detected sample at 
2 = 1.4 — 2.6 in the rest-frame UVJ diagram. The solid line sepa¬ 
rates blue and red star-forming galaxies. Bottom: SFR(Hq:) versus 
M* for the blue and red star-forming galaxies. The log(SFR)- 
log(M*) slope of the blue star-forming sample is steeper than the 
slope of the whole sample, including the red star-forming galaxies. 

galaxies in the MOSDEF sample (Reddy et al. 2015), 
and the CCM curve (Cardelli et al. 1989). The first two 
are stellar attenuation curves while the latter is a line 
of sight extinction curve that we used to correct the Ha 
luminosities throughout this paper. The results are pre¬ 
sented in the bottom panel of Figure 5. The nebular 
color excess (i.e., E{B — F)nebuiar) does not change sig¬ 
nificantly between the three curves, the reason being that 
the shapes of these three curves in the rest-frame optical 
are very similar. The main difference arises from nor¬ 
malization of the curves, affecting the total attenuation 
correction. Reddy et al. (2015) showed that assuming 
the Calzetti curve for the nebular lines as opposed to the 
CCM curve results in higher SFRs, where the difference 
increases with increasing SFR (see Figure 21 in Reddy 
et al. 2015). We see the same trend here in Figure 5, 
namely that the discrepancy between corrected SFRs is 
larger at higher SFRs. Despite the change in SFR, the 
three main sequence slope estimates are consistent with 
each other within their uncertainties. 

5.2. Sample Selection 

Selection effects can have a major effect on the derived 
log(SFR)-log(M,) slope. The samples that are selected 
based on rest-frame UV colors (such as Lyman-break 
selection, Steidel et al. 2003, 2004) are biased towards 
less dusty galaxies. Adding the redder and more dust- 
attenuated star-forming galaxies to these samples results 
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Figure 7. Top: the SFR-M* relation for stacks of galaxies in 
five bins of stellar mass, the difference being in how the stacked 
Ha and H^S luminosities are corrected for the underlying Balmer 
absorption. For illustrative purposes each set of colored points 
is slightly shifted in mass direction. Blue symbols: the Balmer 
absorption in each bin is estimated as the weighted average of the 
individual absorption fluxes. Orange symbols: the same as the blue 
points, but the average is used instead of the weighted average. 
Green symbols: No absorption correction is applied. Bottom: ratio 
of the H/3 absorption flux to the emission flux as a function of 
stellar mass. The absorption flux becomes stronger with stellar 
mass, which results in a more significant correction for galaxies 
with higher stellar masses. 

in shallower slopes (Whitaker et al. 2012b; Speagle et al. 
2014). To test this effect on our sample, we divided the 
star-forming galaxies into a blue bin and red bin, based 
on their rest-frame UV and VJ colors. Previously, we 
used the UVJ diagram to remove the quiescent galaxies 
from our sample (Section 2.2), and now we set another 
criterion on the star-forming part of the diagram to sep¬ 
arate the blue and red star-forming galaxies. The selec¬ 
tion regions are shown in Figure 6. After the removal 21 
galaxies that were classified as red star-forming, the slope 
of the relation increases to 0.75 ± 0.11. These galaxies 
are both more dust-attenuated and older compared to 
the rest of the sample. A steeper slope is expected as the 
red star-forming galaxies dominate the lower right corner 
of the main sequence (those with lower specific SFRs). 
In conclusion, the red star-forming galaxies should be in¬ 
cluded in a complete sample, otherwise the slope of the 
log(SFR)-log(M*) relation may be overestimated. 

Considering that the MOSDFF sample is a near-lR se¬ 
lected sample and because the Ha luminosity is less af¬ 
fected by dust attenuation than the shorter-wavelength 
UV luminosity, this study is less prone to the bias asso¬ 
ciated with UV selection. 

5.3. Balmer Absorption Correction 

Another factor that affects the SFR-M* parameteri¬ 
zation in studies that use the hydrogen Balmer lines is 
the Balmer absorption correction. For the star-forming 
galaxies considered in this study, the difference between 
the absorption corrected SFRs and those that are not 
corrected becomes larger at higher masses because the 
absorption flux, which is produced in the atmosphere of 
primarily A-type stars, increases as a function of stellar 
mass (see the bottom panel in Figure 7). The galaxies 
with higher stellar masses have more mature stellar pop- 
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ulations and hence higher Balmer absorption equivalent 
widths. 

In Figure 7, we show the stacks of SFR(Hq;) in bins of 
stellar mass, corrected for the Balmer absorption in two 
different ways as well as the stacks without any correc¬ 
tions. The absorption correction affects both the value 
of the observed SFR(Ha) and the dust correction (i.e., 
the ratio of Ha to H/3). We used the following correc¬ 
tion method in this study. In each stellar mass bin, the 
Balmer absorption of individual galaxies inferred from 
the stellar population models were averaged, weighted 
by the inverse variance calculated from the errors in the 
individual Balmer emission line luminosities, and added 
to the stacked emission line luminosities. These are blue 
symbols in Figure 7. 

On the other hand, as presented in Figure 7, not cor¬ 
recting for the Balmer absorption (green symbols) results 
in values that are inconsistent with the blue linear fit 
by as large as 9a at log(M*) ^ 10.6. The absorption 
corrections reduce the inferred dust-corrected SFRs by 
values that are typically larger at larger masses. (Here, 
by factors of 1.3, 2.1, and 2.5 at the stellar masses of 
^ 10^°-^, I0^° ®, 10^^-^ Mq.) As a result, the SFRs at 
higher masses are lower once the Balmer absorption cor¬ 
rection is applied. In other words, not correcting for the 
underlying Balmer absorption results in higher SFRs at 
high masses and hence, a steeper slope. 

In conclusion, we emphasize the importance of ac¬ 
counting for the Balmer absorption when using the 
Balmer emission lines to quantify the SFR-M* relation. 
Using the Balmer lines that are not corrected for the 
underlying absorption results in higher SFR values, es¬ 
pecially at the large masses where the Balmer absorp¬ 
tion is stronger, and affects the parameterization of the 
SFR-M* relation. The stacked points in bins of stellar 
mass, for which the weighted averaged Balmer absorp¬ 
tion corrections are applied, are in the best agreement 
(within ~ I — 2cr) with the log(SFR) — log(M*) linear 
regression fitted to the individual galaxies (the blue line 
in Figure 7). 

6 . COMPARISON TO OTHER STUDIES 

The SFR-M* relation has been studied extensively. 
The parameterization of the SFR-M* relation depends on 
the SFR indicator and the typical timescale over which 
it traces star formation. First we consider studies that 
also use Ha to quantify SFRs and therefore provide the 
most direct comparison. Koyama et al. (2013), Sobral 
et al. (2014), and Darvish et al. (2014) used Ha narrow- 
band selected samples at z ~ 2. Such studies, however, 
tend to be biased towards galaxies with larger Ha equiv¬ 
alent widths. These three studies also did not measure 
H/?, and relied on a mass-dependent dust correction from 
Garn & Best (2010). The Garn & Best (2010) relation 
was derived based on a local sample of galaxies but its 
validity has not been verified for z ~ 2 galaxies. The gray 
symbols in Figure 8 are the data points from Sobral et al. 
(2014). At the low mass end, M* < 10^° Mq, the high 
equivalent width bias mentioned above and the effect of 
the improper dust correction are evident. When we ap¬ 
plied an average mass-dependent dust correction based 
on our Balmer decrement measurements to the Sobral 
et al. (2014) observed SFRs, we see a better agreement 
with our results at low masses. The Sobral et al. (2014) 



Figure 8. The comparison of SFR-M* relation for various stud¬ 
ies. The blue points and the blue line are our z = 2.09 — 2.61 
results, using SFRfHo) corrected for dust according to the Balmer 
decrement. The gray squares are from Sobral et al. (2014) narrow- 
band selected star-forming sample, where they used SFR(Ho) dust 
corrected based on the local mass-dependent correction relation of 
Garn & Best (2010). The Whitaker et al. (2014b) and Schreiber 
et al. (2014) relations are based on SFR(UV)-|-SFR(IR). We also 
plotted our results based on the SFR(UV) and SFR(SED) values. 
The relations from the literature are converted to a Chabrier IMF, 
when necessary. 

main sequence is tighter compared to ours, mainly be¬ 
cause the scatter in the dust extinction-M* relation is 
not taken into account. 

Atek et al. (2014) is another Ha study, where they 
used HST/WFC3 near-IR grism spectroscopy for a sam¬ 
ple of emission line galaxies at 0.3 < z < 2.5. They 
used the Ha, H/3, [Oil]A3727, and [Oiii]A5008 nebular 
emission lines as the SFR indicators, dust corrected in 
bins of stellar mass. The log(SFR)-log(M*) slope they 
found was ~ 0.4 at z ^ 2.3, though they also identified a 
Malmquist bias as being a potential issue in their anal¬ 
ysis. The SFR-M* relation was also investigated by Erb 
et al. (2006c) using a sample 114 star-forming galaxies 
at z ~ 2 with Ha spectra corrected for dust extinction 
by the SED E{B — V) and the Galzetti curve. Zahid 
et al. (2012) used the Erb et al. (2006c) sample to fit 
a linear relation to the median SFRs in bins of stellar 
masses. Similar to the other Ha studies, they also found 
a shallow slope of 0.46 ± 0.07. 

We also considered SFR-M* studies where the total 
SFR was estimated as the sum of UV and IR-based SFRs, 
plotting the SFR-M* relations of Whitaker et al. (2014b) 
and Schreiber et al. (2014) on top of our points in Fig¬ 
ure 8. The log(SFR)-log(M*) slopes of these studies are 
generally steeper than our slope. As we used the Balmer 
decrement to correct the observed SFR(Ha) for dust ex¬ 
tinction, as opposed to using the IR data, it is proba¬ 
ble that we missed optically thick star-forming regions, 
particularly in massive galaxies with high specific SFRs. 
This bias may be the reason of our shallower slope com¬ 
pared to the studies with the IR data. There are sev¬ 
eral important differences between the Whitaker et al. 
(2014b) and Schreiber et al. (2014) studies and ours, such 
as the samples, the SFR diagnostics, and the statistical 
methods used for the fit. At this point, it is not trivial 
to reconcile these discrepancies and discuss the agree¬ 
ment or disagreement between these studies and ours. A 
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more informative comparison will become possible once 
we incorporate IR data into our analysis. 

A few studies found indications of a non-linear SFR- 
M* relation (e.g., Schreiber et al. 2014; Whitaker et al. 
2014b; Lee et al. 2015). They suggested that the 
log(SFR)-log(M*) slope flattens at high masses so that 
a single power law cannot explain the SFR-M* relation. 
According to these studies, at higher redshifts the high- 
mass flattening becomes less prominent and at 2 : ~ 2 the 
turning point shifts to M* ~ 10^°'®. A turning point in 
the SFR-M* relation is not obvious from the individual 
detections in our sample, nor from the stacks (Figure 2). 
However, because of the large uncertainty in the most 
massive stacked bin of our galaxies (see the blue and or¬ 
ange points in Figure 7), we can not confidently rule out 
the flattening at the high-mass end. 

7. CONCLUSIONS AND DISCUSSION 

As part of the MOSDEF survey, we used a sample of 
185 star-forming galaxies with Ha and H/3 spectroscopy 
at M* > 10®'® Mq to study the star-forming main se¬ 
quence relation at z ^ 2. The parent catalog is if-band 
(i.e., rest-frame optically) selected and is accompanied by 
ancillary multi-band photometry (Skelton et al. 2014), 
from which the UV-based SFRs were measured. The 
stellar masses were derived by comparing the emission¬ 
line corrected photometric SEDs with stellar population 
models. 

After taking into account the measurement un¬ 
certainties in the SFR and stellar mass, we found 
a measurement-subtracted scatter of 0.31 dex in the 
log(SFR(Ha))-log(M,) relation at 2.09 < z < 2.61, 
which is 0.06 dex larger than what we found based on the 
UV SFRs (see. Table 1). Although in theory the time- 
scale variations in SFR can be traced by using different 
SFR diagnostics, such as UV and Ha, in the SFR-M* re¬ 
lation (Sparre et al. 2015), we argue that in the absence 
of direct measurements of galaxy-to-galaxy variations in 
the attenuation curves and the IMF, the difference in 
the scatter of the log(SFR(Ha))-log(M*) relation and 
the scatter of the log(SFR(UV))-log(M*) relation can not 
be used to constrain the stochasticity of star formation. 
Given these variations, the SFR(UV) scatter could be 
as large or greater than the SFR(Ha) scatter. Theorists 
trying to reproduce the observed scatter in the log(SFR)- 
log(M*) relation should proceed with caution, based on 
the results shown in this paper. 

The scatter in the log(SFR)-log(M*) relation is tighter 
when the SFR(UV) is corrected for extinction by E{B — 
V) derived from the SED modeling or when the SED- 
inferred SFRs are used, because the UV slope regulates 
the color-excess in the SED fitting from which the stellar 
masses were derived. 

In addition to the variations in the dust attenua¬ 
tion, IMF, and metallicity, the scatter in the log(SFR)- 
log(M*) relation may be due to other variables such as 
geometry, size, SFR surface density, and age (e.g., Wuyts 
et al. 2013; Hemmati et al. 2014). The significance of 
the effect of these variables on the scatter will be inves¬ 
tigated in future studies. Although in agreement with 
some simulations (Torrey et al. 2014; Sparre et al. 2015), 
the measurement-subtracted scatter we found is larger 
than some of the other works (e.g., Dutton et al. 2010a; 
Dominguez et al. 2015). The discrepancy may indicate 


that the galaxy accretion and feedback processes are not 
as smooth as predicted by the simulations. 

We found a constant slope of 0.58 ± 0.10 for the 
z = 2.1 — 2.6 sample and 0.65 ± 0.08 when we used 
the whole sample at z = 1.4 — 2.6, for M* > 10®'® Mq. 
Galaxy evolution simulations generally find a close-to- 
unity slope for the log(SFR)-log(M*) relation (e.g., Dut¬ 
ton et al. 2010a; Torrey et al. 2014; Sparre et al. 2015), 
which translates to a stellar mass-independent specific 
SFR (i.e., SFR/M,). Our shallower slope indicates 
higher (lower) SFR at low (high) stellar masses. As 
the slope of the log(SFR)-log(M,) represents the star- 
formation efficiency, one explanation for close-to-unity 
slope in the galaxy simulations is that star formation is 
too inefficient at the low stellar-mass end due to feedback 
processes. The steep slope in the simulations may also 
be caused by a lack of quenching feedback that causes a 
dearth of red star-forming galaxies at the high-mass end 
of the SFR-M, relation. 

The slope of the log(SFR)-log(M*) is influenced by var¬ 
ious observational and measurement effects that could 
explain discrepancies in the slopes reported in different 
observational studies. We demonstrated three of the 
main effects here: dust correction, sample biases, and 
stellar absorption correction (the last one is specific to 
the emission line studies). Gorrecting the Ha emission 
line with the SED inferred E{B — U)sed x 2.27, with the 
assumption of the Galzetti curve, leads to a log(SFR)- 
log(M*) slope of close to unity. A bias against dusty 
galaxies with low specific SFRs, as might be the case for 
UV-selected samples, may also result in steeper slopes 
of the log(SFR)-log(M,) relation. Finally, we showed 
that the Balmer absorption correction has a significant 
effect on the dust-corrected SFRs and not applying the 
correction to the Balmer emission lines could lead to an 
overestimation of the log(SFR)-log(M*) slope. 

The evolution of the slope, normalization, and the 
scatter of the log(SFR)-log(M*) relation is essential to 
our understanding of the processes that govern galaxy 
growth during cosmic time. Once the whole MOSDEF 
sample becomes available, we have access to a large sam¬ 
ple of rest-frame optical spectra for galaxies at 1.4 < 
z < 3.8, which makes it possible to study the SFR-M* 
evolution in redshift, using instantaneous and accurately- 
measured SFRs in a representative sample. Furthermore, 
incorporating IR data into our analysis will allow us for 
a more detailed comparison with the existing studies. 
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